Can dark energy be gravitational waves? 
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The idea that dark energy is gravitational waves may explain its strength and its time-evolution. 
A possible concept is that dark energy is the ensemble of coherent bursts (solitons) of gravitational 
waves originally produced when the first generation of super-massive black holes was formed. These 
solitons get their initial energy as well as keep up their energy density throughout the evolution 
of the universe by stimulating emission from a background, a process which we model by working 
out this energy transfer in a Boltzmann equation approach. New Planck data suggest that dark 
energy has increased in strength over cosmic time, supporting the concept here. The transit of these 
gravitational wave solitons may be detectable. Key tests include pulsar timing, clock jitter and the 
radio background. 



I. INTRODUCTION 

Dark energy was originally detected as accelerated ex- 
pansion seen in the distance scale for supernovae of type 
la (Schmidt et al., 1998, Ricss et al. 1999, Perlmutter et 
al. 1999; for a review see Frieman et al. 2008). Many 
suggestions have been made about what dark energy is, 
what its strength is, what its time evolution is, and what 
possible further observational results are. 

The idea that dark energy is gravitational waves may 
explain its strength and its time-evolution. One possi- 
ble concept is that dark energy is the ensemble of co- 
herent bursts (solitons) of gravitational waves originally 
produced when the first generation of super-massive back 
holes was formed (Caramete & Biermann 2010); the en- 
ergy density of such solitons would suffice within the un- 
certainties. These solitons get their initial energy as well 
as keep up their energy density throughout the evolu- 
tion of the universe by stimulating emission from a back- 
ground (Biermann & Harms 2013). Our model of the 
background metric resembles the Randall-Sundrum ideas 
(1999a, b) but is time-dependent, and describes the en- 
ergy flow from the background (strong-gravity) brane 
to our world (weak-gravity) brane. Planck data sug- 
gest that dark energy has increased in strength over cos- 
mic time (Planck 2013 XVI), as predicted by our model. 
Gravitational waves were far below today's dark energy 
at the epoch of early nucleosynthesis and of the forma- 
tion of the microwave background ripples (as summarized 
in Ligo+Virgo-Coll. 2009), both much earlier than the 
likely formation epoch of the first generation of super- 



massive black holes. Our model is also consistent with 
early star formation, as we argue below. The transit 
of the gravitational wave solitons postulated here may 
be detectable. We discuss the predictions briefly below 
and elsewhere. We focus on the Boltzmann equation ap- 
proach, working out the energy transfer from the strong 
gravity background in stimulated emission. 



A. Basic ideas 

Inspired by Bekenstein's (1973) considerations we 
posit: When the first generation of super-massive black 
holes was formed, each produced a coherent burst of 
soliton-like gravitational waves which combine to give a 
total energy of order 
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In the following we also call this an ensemble of soli- 
ton waves, or shell fronts. Nbh,o is the original comov- 
ing density of super-massive black holes. Today super- 
massive black holes have a density of io _1 - 7±0 - 4 Mpc -3 
above Mbh = 3 ■ 10 6 M Q (Caramete & Biermann 2010); 
assuming that they grow by merging, and allowing for 
statistical and systematic errors, an original comoving 
density of Nbh.o = 1 Mpc -3 seems possible. This co- 
moving density is the density black holes had at the 
beginning, so transposed to today without change in 
their numbers per comoving volume. The data sug- 
gest that there was a generation of first super-massive 
black holes with a mass between Mbh ~ 10 6 Mq and 
Mbh ~ 10 7 M Q . The original black hole mass may be 
~ 3 ■ 10 6 Mq considering (i) the black hole mass function 
(Greene et al. 2006, Caramete & Biermann 2010), (ii) the 
instability of massive stars (Appenzeller & Fricke 1972a, 
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b) in an agglomeration picture (Spitzcr 1969, Sanders 
1970), and (iii) the observed black hole in our Galactic 
Center (e.g. Eckart et al. 2005). The redshift of creation 
z+ may be large, as formation of massive stars may begin 
at redshift 80 (Biermann & Kuscnko 2006). Redshifts 

from about 30 to 50 allow a quantitative interpreta- 
tion of the data of dark energy. At the original density of 
black holes adopted here redshift 50 is consistent with the 
mass of ~ 3- 10 6 M Q , and redshift 30 would imply 

Mbh ~ 10 7 Af Q , in either case to make the estimate 
consistent with dark energy today. 

What is the motivation for considering gravitational 
waves? Bekenstein (1973) wrote about the entropy of 
the universe: "... we must regard black hole entropy 
as a genuine contribution to the entropy content of the 
universe" . However, entropy is also information, and in- 
formation must have a carrier. A natural suggestion is 
that this carrier is gravitational waves, with an energy 
commensurate with the black hole scale. This suggestion 
is consistent with the fact that cosmological black holes 
are not in thermodynamic equilibrium, and therefore the 
entropy associated with such black holes should be de- 
scribed by statistical mechanics as advocated in Harms 
and Leblanc (1992, 1993). This speculation immediately 
gives 
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where Ngw,q is the number of gravitons at the forma- 
tion of the black hole. For M BH = 3 • 10 6 M Q this is 
Ngw,o — 10 90 . Mbh is the original mass of the black 
hole, vn pi is the Planck mass, c is the speed of light, and 
Gjv is Newton's constant of gravity. Eqw is the average 
graviton energy given by 
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4D strong gravity brane and our 4D world weak-gravity 
brane. 

ds 2 = - e (V0 m */* c 2 dt 2 + eitY W du 2 

+e (i-Kf ) n )a*/«te e -(T) fc ( 1 -D dxt dx* , (5) 

where i = 1,2, 3, Th is the Hubble time, I = Ipi is 
the Planck length, Tpi = lpi/c is the Planck time, u is 
the coordinate in the fifth dimension, and the remaining, 
non-coordinate quantities are arbitrary parameters. Al- 
though the five-dimensional covariant divergence of the 
energy-momentum tensor does not vanish everywhere, it 
vanishes on the weak-gravity (weak) brane (it = 0), and 
it is approximately zero on the strong-gravity (strong) 
brane (u = I) for very small dimensionless ratios /3/th 
and ip/rn with ip < and /3 > 0, with the conditions a) 
r H {ip + » 1, b) 2(6 - l)/a > t h I4> - 1, and 

c) a 2 = 3. goo in the metric above then defines the con- 
fining potential for the strong brane. Gravitons from our 
brane stimulate transitions between bound states on the 
strong brane, resulting in the emission of gravitons onto 
our brane. This metric describes a weak brane for our 
world which is expanding with time, and a strong brane 
which is contracting with time, albeit very slowly for the 
latter brane. The 5D-cosmological constant measured on 

the weak brane is £l wea k = — \ Tir ) ' 

In our model the gravitons at high energy in the back- 
ground obey a Planck-like distribution with Planck tem- 
perature. We adopt the point of view that the Planck 
scales are limits: nothing can go below Planck time 
and Planck length, and no single particle can go beyond 
Planck energy, in any frame. The strong brane is stable 
against collapse, since given a Planck spectrum for any 
wavelength A the free-fall time scale 77 / is always either 
equal or longer than the pressure wave time scale r s . 

r ff = r «fe) 3/2 ^ r «fe) = r - 



This gives a graviton energy of Eqw — 10 30 erg for this 
black hole mass. The entire energy content then is 
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We picture this as a coherent burst of gravitational 
waves, or a soliton wave, ejected at formation of the black 
hole. It is clear from the considerations above that we 
are not using the weak-field approximation. Multiply- 
ing with the original density of super-massive black holes 
reproduces our estimate above. 



II. ENERGY TRANSFER BY STIMULATED 
EMISSION 

In the following we use a particle-wave duality for 
gravitons at high energy, which thus associates a 
wavenumber k/h and a corresponding length-scale A to 
each spatial direction, and we assume localization is pos- 
sible to about a wavelength. 



A. Direct derivation from Boltzmann equation 



B. Background model 

In our model for the background, which has some simi- 
larity to the Randall- Sundrum (1999a, b) ideas, we iden- 
tify a possible local metric to describe a 5D world with a 



The distribution function Af(k, t) satisfies the equation 

I f ri 3 k' f cd 3 p r cd 3 p' r d 3 k" 
kJ (27r) a 2fe' J (2tt) 3 2 E{p) J (2x) ;1 2 E(p') J (27r) :i 2 k" 
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|M| 2 2 n 4 S 4 (K + Q-K'- K" - Q') 2 tt 3 5 3 (k< - k 1 ') 
(Af(k', t)M{k", t) g b (p', t) (1 + Af(k, t)) 

-Af(k,t)g b (p,t)(l+N(k\t))(l+Af(k",t))) . (6) 

where K = (k,k), Q = (E,p), and Af(k,t) is the dis- 
tribution of allowed states for gravitons with momenta 
k = \k\ and p = \p\ at time t on the shell. The (^-function, 
S 3 (k' — k"), has been inserted to impose coherence of the 
outgoing gravitons. \M\ 2 is the matrix element squared 
for the quadrupole emission of a graviton of 4-momentum 
K" . gb{p, t) is the occupation number distribution of the 
background particle sea. R(t) = (1 + + z) is the 

scale factor for an expanding universe. 

The Boltzmann equation for Af(k, t) to lowest order in 
the expansion of the 4-dimcnsional (5-function is 
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where the factor k is given by 
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and \M\ 2 is as stated above. 

The product N{k, t) {Af{k, t) + 1) can be expressed as 

- N{k, t) (Af(k, t) + l) = ^T g t) , (9) 

where T g = T gQ / R{t) and k is the magnitude of the gravi- 
ton momentum. This replacement is valid for a Planck- 
like distribution on the shell. After making the replace- 
ment k = k/R(t), the equation for Af(k,t) becomes 
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or in terms of the frequency of the wave at emission isq 
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The solution of this equation can be written in terms of 
dimensionless parameters x and y, where X = ^~jr , and 



y = Jl ^r^dt' a 
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arbitrarily set y = and will discuss the possible signifi- 
cance of the term elsewhere. 

The rate at which energy is created can be calculated 
from the expression for J\f in cqfTT] The rate of energy 
creation per one original graviton (of the total number 
Ngw.q) is 

d < E > 



dt 



The total rate of energy creation is 
d < E T > 



dt 



N GWfi R(t) 4 Ko cA 
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A = J$% J 'vf-^dv,, and k B T g0 = 

In analogy with electrodynamics we adopt the point of 
view that the matrix element squared for the emission of 
a quadrupole wave from the background evolves as 



(15) 



This dependence on redshift can be thought of as one 
factor of (1 + z) as arising from the redshift dependence 
of energy over mass c 2 , and the other three powers as 
(1 + z) 3 cx 1/R 3 arising from the spatial dependence of 
quadruple emission. The overall time scale can reason- 
ably be expected to be H(z)~ 1 . We expect A to be of 
order unity. 

This allows the expression for rf< ^ T> to be consis- 
tent with the observed energy density and the interaction 
term kq to evolve with red-shift as 
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Inserting this parameter dependence then leads back, 
to within a factor of order unity, to the result we were 
seeking, 



(17) 



After integrating we obtain with this redshift dependence 
a constant dark energy density as in eq. [TJby multiplying 
by the redshift evolution of black holes Nbh,o (1 + z) 3 . 



with X = ^/{x 2 + 2y}. When y reaches the same 
order of magnitude as x 2 , the number density J\f(x, y) be- 
gins to deviate from the Planck-like distribution; this is 
to be expected as soon as the energy becomes low enough 
for mass to appear. This leads to the speculation that 
this term might correspond to the Higgs mass. Here we 



B. Heuristic derivation 

The microscopic picture derives from the resonance ar- 
gument in momentum space between our weak gravity 
brane and the strong gravity brane. In analogy to the 
terms in the Kompaneets (Kompaneets 1957, Zeldovich 
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1975) equation describing stimulated emission for pho- 
tons, derived from the Boltzmann equation, there are 
four steps necessary: We need to 1) define the reference 
frame for the stimulation; 2) define the phase space den- 
sity, which then occurs squared in the relationship for 
stimulated emission; 3) determine the number of new 
gravitons created; and 4) fix the characteristic energy 
and rate. 

In the determination of a possible frame of reference 
the key point is that the coherent wave can be localized 
only to a finite length. We start with the wavelength 



of the resonant wave: \qw — 3 Ipi 



Mr 



For a very 



large number of gravitons Nqw, however, we obtain a 
precision in localization of 
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A localization of a signal with uncertainty implies a ve- 
locity < c. Nominally, a factor of (1 + z)/(l + Z+) < 1 
ought to be there, but as in our approach no length 
scale can be shorter than the Planck length is any frame, 
the Planck length itself is the limit. A one-step ran- 
dom walk in space analysis to describe the uncertainty 

yields a Lorentz factor given by ^ _ 1 
r{z,z*) = I 

to observer. Therefore T = r[z,Z+)/lpi. Space un 



- 2r^, where 
is the distance scale from source 



H(z') 

Therefore T = r(z,Z ic )/lpi 
certainty implies going forward or backward, time un- 
certainty implies going forward or not going forward. 
For the time uncertainty we obtain an expression for the 
Lorentz factor with a square-root. We use H(z)~ 1 as an 
approximation for the time rather than its redshift inte- 
gral, since it allows a simpler final expression, and is a 
good approximation for times near today. Combining two 
Lorentz factors allows either the addition or the subtrac- 
tion of the velocities, and here only the subtraction gives 
reasonable numbers - unreasonable is an implied energy 
of a single graviton beyond the Planck scale. Therefore 
using both space and time gives finally the Lorentz factor 
of a shell front 



2 Ipi 



(19) 



For a high original redshift z+ looked at from today (z = 
0) this Lorentz factor is very large, typically of order 10 31 . 

In momentum space the gravitons are coherent within 
a coherence patch, a fact we will use below. A coherence 
patch is defined by the wavelength of the gravitational 
wave or graviton in the two directions perpendicular to 
the direction of the highest momentum component. In 
momentum space the gravitons are highly anisotropic, 
but we are concerned with the rate of interaction with the 
background, and thus have to multiply the cross-section 

with the Lorcntz-factor T to obtain (^Ngw ~§^p2 > 

where we have normalized this surface density to the 
Planck area. 



We can derive this critical phase space density by not- 
ing, that multiplying this term with the coherent area 
7rA 2 , the collision velocity in the comoving frame of c, the 
Planck time Tpi, the density scaling of the background 
Ipf, the interaction probability V, and the fraction of res- 
onant phase space in the Planck distribution of the back- 
ground (mpi/Mgjj) 3 r/(r mpi/Afsg) gives unity for ex- 
actly V = 1, so stimulating the energy transfer of one 
graviton to our weak brane. This is consistent with 
Fermi's Golden rule in which the rate of energy trans- 
fer depends on the number of states linearly. 

In a numerical example of a black hole mass of 3 • 
10 6 Mq, a creation redshift for the first generation of 
z* = 50, and a scale R s equal to the luminosity distance 



<*pi 

4tt_R2 



cIl the critical phase space density \Ngw 

of order unity. We note that the luminosity distance is 
in our general case = r(z,z*)(l + z*)/(l + z), just 
taken from redshift z to redshift z*. 

The number of coherence patches is, again in the co- 
moving frame, ^ , which is just the total area 

divided by the lateral area of a coherence patch. Again 
we note that we are talking about a rate, and so get the 
extra factor of V in the shell frame. 

We use stimulated emission, so the energy of the single 
new graviton En pulled from the background has to be 
in resonance: Er = EqwT. The rate is just one 
graviton per Planck time in the proper frame (rate) = 

Tp^ . These expressions are again in the shell frame. 

Out of these four terms we can form an expression 
for stimulated emission. This equation for the stimu- 
lated emission is constructed in a familiar non-linear form 
(Kompaneets 1957) for stimulated emission and can be 
written in the form of energy E per unit time created in 
the form of resonant gravitons (Bicrmann & Harms 2012, 
2013), putting all the terms together: 



dE - ( N apl ^ 
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in the comoving frame with proper time r on a spherical 
surface of radius R s with graviton number Ngwi now 
including its redshift dependence. In eqUHlwe have dif- 
ferentiated between the Lorentz factor for the surface, 
and the energy transformation, since in a more general 
application they may not be the same. Te is the Lorentz 
factor as applied to the transformation of energy. Eqw is 
the energy of the gravitons on the surface in the observer 
frame, and A ~ ch/Eow is the corresponding length 
scale. Td is a Lorentz factor applied to the surface den- 
sity relative to the observer and to the lateral surface 
density in terms of coherence patches of area ir A 2 . In the 
following analysis we identify these two Lorentz factors 
with the Lorentz factor T derived above. 

As a check of this equation we realize that the number 
of coherent patches is to within a factor of the Lorentz 
factor r again the critical phase space density, but in- 
verted. This implies that we can simplify the equation 
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above and thus resurrect the interaction rate including 
the energy of the resonant graviton produced in stimu- 
lated emission. 



C. Consistency 

We can take the time derivative of the existing dark 
energy and compare with the derivative of the equation 
derived above heuristically for the stimulated emission: 
4:(NgwEgw) f° r which we obtain a solution with the 
redshift dependency Nqw OC ( ) ■ Using E@yy OC 

^ we obtain § M BH C 2 (^) 3 H(z). 

Comparing cq|17lwith the stimulation expression eg 1201 
gives a mismatch by (tt/3) 1 ^ 2 ~ 1.02, which could easily 
be incorporated into the definition of the critical phase 
space density for instance. So the key result of this paper 
is the redshift dependence for New, which gives a quan- 
titative solution to the non-linear stimulation expression 
cql2"01 The initial level of Nqy/ is Bekenstein's (1973) 
entropy S/Hb- 

D. Caveats 

First, we have used the assumption that there is actu- 
ally a single identifiable first generation of super-massive 
stars, which evolve into black holes, all of the same mass. 
However, as can be easily ascertained, any later deposi- 
tion of gravitational wave energy is weakened relative to 
the first by the product of the density of those black holes 
Nbh,o and the redshift factor (1 + z+) 3 . An alternate 
possible picture would be the formation of degenerate 
dark matter stars, which get eaten up from the inside by 
stellar mass black holes (Munyaneza & Biermann 2005, 
2006) to produce super-massive black holes. A mass of 
order 3- 10 6 M Q (Wang & Biermann 1998) allows the big 
mass black holes to form relatively early in the universe 
as observed; any original mass significantly less would 
not allow sufficient time. 

Second, we have assumed, that black holes reach their 
high masses observed more by merging than by accre- 
tion, since with such a model we obtain the relatively 
high early black hole densities necessary. The electrody- 
namic emission of accreting black holes integrated over 
time tends to exceed what is implied by cosmological sim- 
ulations (Caramete & Biermann (2010)). All the energy 
deposition "visible" in the sum is accounted for by large 
scale structure formation. 

A third assumption is that stellar black holes do not 
participate significantly. The summed contribution of 
stellar mass black holes might be expected to exceed the 
contribution from a single central super-massive black 
hole. However, stellar black holes do not provide a criti- 
cal phase space density for their ejected shell of gravitons, 
and so "drop out" . This cut-off is clear from the deriva- 
tion. 



III. PREDICTIONS 

First, dark energy should be detectable by pulsar tim- 
ing (e.g. Kramer 2010). Our proposal says that dark 
energy arises from coherent bursts of gravitational waves 
produced by forming the first generation of super-massive 
black holes. This implies that at the properly redshiftcd 
frequency of these black holes the gravitational wave 
background should reach a value close to unity relative 
to critical density, with a Planck spectrum. 

Second, the passing of the soliton-like gravitational 
wave shells should be detectable as jitter in clocks (e.g. 
Predehl et al. 2012), given a substantial improvement 
over current observations. 

Third, the original creation of the super-massive black 
holes should be detectable through their extended rem- 
nants in the radio, X-ray, 7-ray and neutrino background, 
and in ionized molecular hydrogen absorption. It may 
have been detected already in the radio (Fixsen et al., 
Kogut et al., Seiffert et al. 2011; Condon et al. 2012). 

Fourth, the fragmented gaseous shell around the rem- 
nants could be the birth of bulge-less disk galaxies (Ko- 
rmendy et al. 2010, 2011, Kormendy & Bender 2011). 
Such a population of huge shells may be consistent with 
the density and masses of such galaxies. 

Fifth, intergalactic magnetic fields may be initialized 
at that epoch as well. Massive stars rotating due to tidal 
forces form magnetic fields (e.g., Biermann & Schliiter 
1951), and turbulent motions with rotation will enhance 
these magnetic fields by the dynamo effect. 

IV. CONCLUSIONS 

In this paper we have presented a simple analytic 
picture for the origin of dark energy, based on the 
idea that in their formation super-massive black holes 
eject a coherent burst of gravitational waves, which get 
about (1/2)Mbhc 2 from the background by stimulated 
emission. The acceleration equation now takes the en- 
ergy transfer ((4ttGn)/3)(^pde) of dark energy into ac- 
count with (d 2 R/dt 2 ) = -((47rGjv)/3)(p + 3P/c 2 ) + 
{{AttGn) /3){Apde)- The factor of 4 in the energy trans- 
fer term derives from the compensation of the energy loss 
of a relativistic gas running with (1 + z) A . Interpreting in 
the limit p as dark energy, and using the e.o.s. of gravi- 
tational waves of P = (l/3)pc 2 we resurrect the normal 
equation with an apparent negative pressure, mimicking 
Pde = —PDEC 2 - The background is conceived as simi- 
lar to the model by Randall & Sundrum, but variable in 
time. 

Some tests of these ideas are i) pulsar timing to de- 
tect the gravitational wave background which constitutes 
dark energy; ii) a jitter in clock timing associated with 
the many coherent bursts of gravitational waves going 
past; and, iii) a radio background associated with the 
remnants of the formation of the first generation of black 
holes. 
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